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ABSTRACT 

The evolution of small scale cosmological perturbations is carefully re-examined. 
Through the interaction with photons via electrons, baryon perturbations show inter- 
esting behavior in some physical scales. Characteristic features of the evolution of 
baryon density fluctuations are discussed. In CDM models, it is found a power-law 
growing phase of the small-scale baryon density fluctuations, which is characterized 
by the terminal velocity, after the diffusion (Silk) damping and before the decoupling 
epoch. Then, a transfer function for total matter density fluctuations is studied by 
taking into account those physical processes. An analytic transfer function is presented, 
which is applicable for the entire range up to a solar mass scale in the high— z universe, 
and it is suitable also to the high baryon fraction models. 
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1. Introduction 

The structure formation in the high-z universe is one of the most important problem in the 
fields of cosmology and astrophysics. The structure formation of the Cold Dark Matter (CDM) 
dominated models is well motivated from the recent cosmological observations, e.g., from the mi- 
crowave background anisotropics and the large scale structure of galaxies ( [White, Scott, &: Silk 



1991 podelson. Gates, fc Turner 1996| ; [Peacock fc Dodds 1995| ), although it may require some 



modifications, i.e., inclusion of cosmological constant, and the tilted initial power spectrum with 
the gravity wave mode (White, Scott, Silk, & Davis 1995). We expect that the power spectrum 
of density fluctuations on scales ^ IMpc will be measured precisely in near future, e.g., from the 



2DF survey, the Slon Digital Sky Survey (SDSS) project (see, e.g., Strauss 1996; Loveday 1996) 



and the future cosmic microwave background experiments with satellites ^ (see, e.g., Zaldarriaga, 



Spergel, & Seljak 1997; Jungman, Kamionkowski, Kosowsky, &: Spergel 1996^ . These observations 



will severely constrain the cosmological parameters and the initial density power spectrum. 

On the other hand, the inflation scenario provides a successful mechanism to explain the origin 
of cosmological density fluctuations. It predicts almost scale invariant and Gaussian fluctuations 
( Bardeen, Steinhardt, &: Turner 1983[ ). In this scenario, we expect that the initial density power 



spectrum extends up to very small scales. The evolution of density perturbations due to the 
fluctuations on very small scales is especially important for the early formation of the bound objects 
such as population III stars or primordial sub-galaxies for the hierarchical clustering scenario of 
the structure formation. By studying this evolution up to the non-linear regime with taking into 
account the heating and cooling processes, we would be able to learn the formation epoch and the 



formation process of the initial cosmic objects in the high-z universe (e.g., Ostriker & Gnedin 1996 
Gnedin k Ostriker 199^ ; Haiman k Loeb 1997). 



In this paper, we carefully re-examine the evolution of the small scale cosmological perturba- 
tions in the linear regime as an initial condition of the structure formation. In cosmological models 
with the high baryon fraction, the evolution of baryon density fluctuations affects that of the total 
matter fluctuations. This derives an alternation of the standard matter transfer function of the 
total matter density fluctuations. In particular, the power spectrum shows the damping on scales 
which are smaller than the Jeans scale at the decoupling epoch. This damping is caused by the 
acoustic oscillation of baryon perturbations. After describing these characteristic features of the 
evolution of baryon density fluctuations, we re-investigate the transfer function for the total matter 
density fluctuations. 

In § 2, we give a brief review of the equations which describe the evolution of the cosmological 
perturbations. In § 3, we describe various physical scales which are relevant to the characteristic 
features of the evolution of baryon density fluctuations. The evolution of baryon density fluctuations 



^ Planck home page 
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before and after the decoupling are discussed in § 3 and § 4, respectively. The total matter transfer 
function is examined in § 5. In § 6, the accuracy of the transfer function is shown by calculating 
statistical quantities such as Cg. A physical implication of our results is also demonstrated by 
computing the formation rate of small bounded objects in the high-z universe based on the Press- 
Schechter theory. § 7 is devoted to summary and discussions. In Appendix A, we summarize the 
previous results on the CDM perturbations in small scales by Hu k, Sugiyama (1996) (hereafter 
HS96). In Appendix B, we investigate the effect of the baryon thermal pressure on the transfer 
function on small scales after the decoupling. 

Throughout this paper, we work in units where c = h = ks = And we assume Tq = 2.726 K 
as the cosmic microwave radiation temperature at present unless we explicitly include this valuable 
in equations. 



2. Review of Formalism 

Let us first summarize the equations which describe the evolution of the cosmological density 
perturbations on small scales. We write the perturbed metric in the Newtonian gauge, 

ds2= (^^y(-(l + 2*)dr/2 + (l + 2$)dx2), (1) 

where ^ is the perturbed gravitational potential, $ is the curvature perturbation, a is the scale 
factor and the suffix indicates the present value. As we are interested in the small scale cosmo- 
logical perturbations, we can assume the geometry of the universe to be flat. We take into account 
the effect from the curvature term and the cosmological constant of the universe only near the 
present epoch through the change of the expansion rate. Then, in the high redshift universe, the 
Priedmann equation is written as, 

= (^^y = (^y^^^n,Hi (2) 



a a J \ a J Oeq + ag 

where the dot denotes ?7-differentiation, Hq is the Hubble constant, and Oeq is the scale factor 
at the matter-radiation equal-time. The Priedmann equation can be integrated as keqri/2\/2 = 
1 -I- a/tteq — 1), where we define keq^ = 2fio-f^o ('^o/'^eq)) ^'^d 



fceq = 0.0948 X - fMoh^Mpc-\ (3) 

— = 4.04 X 10^(1 - ^)0o/l^ (4) 

with being the neutrino fraction in the energy density of radiations. Por the massless standard 
neutrino model with three species, fi, ~ 0.405. 



The evolution of photon and neutrino radiation fields is described by the the Boltzmann 
equation. The linearized evolution equation of the k-th Pourier mode of the photon temperature 
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perturbations is (see e.g., Kodama &: Sasaki 1986; Hu & Sugiyama 1995, hereafter HS95), 

G + ikfiie + ^r) = -6 + f (eo - e - ^G2P2(/u) - ifiV^) , (5) 

where /i = 7- k/|A;|,f = UeCrTa/ao, Ue is the free electron number density, ctt is the Thomson cross 
section, Pi^n) is the Legendre function, Vb is the velocity of baryon perturbations, 0; is defined by 
the multipole moment expansion of 6(77,^,7) = S/^o ^K^?; ^) ^^'^ 7 is the directional 

vector of the photon momentum. We can solve the equation using the multipole expansion method. 
The perturbations for the neutrino radiation fields obey to a similar equation. Effects of neutrino 
perturbations, however, are not so important for the structure formation on small scales. It should 
be taken into account only around the epoch of the horizon crossing (HS96). 

In the fluid approximation, the evolution of CDM density fluctuations is described as 

5c = -kVc - 3$, (6) 

Vc = --Vc + k^, (7) 
a 

where 6c and Vc are the fc-th Fourier mode of the perturbations of the CDM density contrast and 
the velocity, respectively. 

In our gauge, the perturbed Einstein equation gives the extended Poisson equations, 

2 



\aoJ 

and 



= 47rG(^^^ Pt(5t + 3^(1 + u;t)A;"Vt) , (8) 
(-) ptUt, (9) 



2 

/c2($ + ^) = -8ttG{ 

where pril + wt)Vt = T,xiPx wt = Pt/pt, Pt^t = T,xPxSx, the stress anisotropy is 

IIj' = 12(02 + -^2)75, 02 and A'^2 are the quadrupole anisotropies of the photon and the neutrino 
temperature perturbations. 

We next consider the baryon-electron system. Electron(e), neutral and ionized hydrogen(H) 
and helium(He) are the particle species of the baryon-electron system. The number densities for 
each species are written as 

Tie = Xe(l - y)nb, nH=(l-yp)nb, TlHe = ^"b, (10) 

where yp is the primordial helium mass fraction, tih and nue are the number densities of hydrogen 
and helium, respectively, = nH+4nHe is the total baryon number density, and Xe = rif,/ (ren+'^He) 
is the electron ionization fraction. We take a single fluid approximation for this baryon and electron 
system. Then the evolution equations of perturbations of the baryon-electron fluid can be written 
as 

6b = -kVb - 3$, (11) 

Vb = --Vb + kc,^6b + k^ + f (01 - Vb)/R, (12) 
a 
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where Cg is the sound velocity of the baryon-electron fluid, and R = 3pfe/4/9^ = 3r2b/40o(l — 

In the previous paper (Yamamoto, Sugiyama, & Sato 1997, hereafter YSS), we gave useful for- 
mulation for the small scale baryon perturbations. In particular we investigated how the Compton 
interaction between the electrons and the background photons determines the sound velocity of 
the baryon perturbations after the recombination. According to the previous result, the effective 
sound velocity of the baryon-electron fluid after the recombination can be determined by solving 
the dispersion relation, 

- i-qELoiu;^ - Cf^/c^) + uj'^ - Ce^k^ = 0. (13) 

Here the real part of uj/k can be regarded as the effective sound velocity, i.e., Cs = Re[a;/A:], t/e is 
the energy transfer time of the Compton interaction, 

8 a Xe(l -yp/2)cJTP7 

= 7 ^ ^ ' (14) 

3 «0 me(l + Xe- (Xe + 3/2)yp/2j 

Cf^ = 5Po/3/56, Ce^ = Po/pb, Pq ^nd pb are the pressure and the energy density of the baryon-electron 
fluid, respectively. Note that Cf is the sound speed for an adiabatic process and Cc for an isothermal 
process. The effective adiabatic index 7 can be defined as 7 = jc^ . As is clear from equation 
( [I^ ) , the effective sound velocity Cg (adiabatic index 7) is ruled by the ratio of the Compton energy 
transfer time scale t/e to the sound oscillation time scale. Namely, if t^eoj <C 1, then lo c^k and 
7 ~ 1. On the other hand, if ^^uj ^ 1, w ~ Cf A; and 7 ~ 5/3. Therefore the sound velocity depends 
on the wavelength of the perturbations, and 7 is changed from 1 to 5/3 as the universe expands. 



3. Evolution of the baryon perturbations before the decoupling 

In this section, we describe the evolution of the cosmological baryon perturbations on small 
scales before the decoupling. The decoupling epoch (or the drag epoch) r/^ of baryon perturbations 
is defined by Tdirfdj^o) = !> where 

Td{Vi,r]2)= dr]f{rj)/R{i]) . (15) 

One can find the analytic fitting formula of this epoch in HS96. Note that this epoch is not 
necessarily equal to the last scattering epoch of photons which is defined by T{rirec,'>lo) = 1- 

From the Compton interaction between electrons and photons, there are interesting varieties 



of physical scales for the baryon-electron fluid. In the previous paper ( YSS ), we summarized the 
details of the physical scales which are relevant to the evolution of baryon fluctuations. In this paper 
we focus on the small scale density perturbations, and describe the evolution using the physical 
scales. We define kP^y^" {XP^y'' = iTr/k^^y'') physical wave number (wavelength), and k (A) as a 
comoving wave number (wavelength). They are related as A; = (a/ao)/cP^^'^ { ^ = {ao/a)XP^y^ ). 
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3.1. Horizon crossing and Jeans oscillation 

The first characteristic scale is the horizon crossing scale in the evolution of perturbations. 
We define the horizon crossing scale by l/fcn = and the corresponding baryon mass scale by 
Mh = ATTpb{i:/k^^y/?,. Then we have Mr = 7.9 x 1030(1 + z)-3(l - f^)^/^Q^h'^ Mq, in the 
radiation dominated stage. According to this definition, the horizon crossing occurs at the redshift 

1 + ZH = 2.0 X 10lO(l - /,)l/2(17b/l2)V3(^/^^)-l/3 ^ (ig) 

for perturbations with the baryon mass scale M. Since we have assumed that the perturbations 
cross the horizon in the radiation dominated stage, this expression is applicable for M <C 1.2 x 
1017(1 - /,)-3/2j7b/i2(J7o/i2)-3 Mq. 

The next important scale is the Jeans scale. Since the coupling between baryons and photons 
is tight due to the Compton interaction sufficiently before recombination, the photon-pressure 
prevents perturbations from collapsing. Therefore the baryon-photon perturbations oscillate as 
an acoustic wave inside the Jeans scale. We define the Jeans wavelength (wave number) before 
recombination by Aj*^-^*^ = 2T:/k^^^ = (7rCs^/G(pm + P'y)Y^'^ where Pm = Pb + Pc and pc is the 
energy density of CDM. Defining the Jeans mass as Mj = 47rpf,(Aj'^^'^/2)'^/3 = 47rp5(7r/A;j'^^^)'^/3, 
we get 

Mj = 8.6 X 10^^(1 + z)-^Qy,h'^ Mq, (17) 
in the radiation dominated stage. Then the Jeans crossing occurs at redshift, 

1 + 2j = 9.4 X 10^{Q],hY^^{M/MQy^^^ , (18) 

for the perturbations with the baryon mass M, which is right after the horizon crossing epoch. 
Since we have assumed perturbations cross the Jeans scale in the radiation dominated stage, this 
expression is applicable for M <C W^^Q^h'^ i^oh'^)~^ Mq. When the coupling between baryons and 
photons is tight, the baryon-photon system behaves as a single fluid. The oscillation is expressed 
in the analytic form, 

'^b = ^(l + ^/.) 'f(0,fc)(l + i?)-i/4cos(A;rs), (19) 

for adiabatic perturbations (HS95;HS96), where rg = /J' di]' /y/3{l + R), and $(0,A;) is the initial 
curvature perturbation. 

Figure |l] shows the numerical evolution of baryon and CDM density fluctuations with various 
baryon mass scales, IO^^^Mq, 10^ Mq, IO^Mq, and 10^ Mq. It is apparent from the figure that 
baryon density fluctuations start oscillating after the Jeans crossing. It is also well known that CDM 
density fluctuations grow logarithmically after the horizon crossing in the radiation dominant regime 
(see e.g., HS96). Moreover, this figure shows some prominent features of the evolution of baryon 
density fluctuations, i.e., damping, re-growth and catching up with CDM density fluctuations. 
These are discussed in the following subsections and §4. 
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Fig. 1. — Typical numerical evolution of baryon and the CDM density fluctuations, (solid lines) and Sc (dashed 
lines), with baryon mass scales IO^^Mq, IO^Mq, W^Mq, and W^Mq, respectively. These mass scales correspond 
to the wave numbers k — 9.65 x 10~^, xlO", xlO^, and xlO^Mpc"^, respectively. Note that here we employ the 
gauge to set the total matter center-of-mass frame, instead of the Newtonian gauge for the convenience of numerical 
calculations. The behavior of fluctuations inside the horizon does not depend on the gauge choice. We chose the 
cosmological parameters, h — 0.5, f2o ~ 1.0, and f2b = 0.1. The dotted line shows the decoupling time of baryons and 
photons. We see the power-law growing phase of baryon density fluctuations after the diffusion damping before the 
decoupling. After the decoupling the baryon fluctuations catch up with CDM fluctuations through the gravitational 
free fall, unless the perturbations are smaller than the Jeans scale (see panel (d)). 

3.2. Diffusion damping and the breaking scale of the tight coupHng approximation 

It is well known that the photon diffusion process erases the baryon-photon perturbations 
on small scales (Silk 1968; Sato 1971; Weinberg 1971). This damping length is the random walk 
distance of a photon scattered by electrons through the Compton interaction. The mean free path of 
a photon is 1/f. The integrated distance that the photon proceeds is the horizon scale rj in terms of 
the comoving length. Hence this damping scale is roughly expressed as / ~ ^frifr. In a more detailed 
treatment, we define the damping scale by k^, where fc^^ = /J' driiP? + 8(1 + i?)/9)/(6f (1 + Rf) 
(see e.g., HS95). It is known that the exponential damping factor, exp[— fc^/Zc^,] is multiplied to 
the right hand side of equation (|l^) to describe this damping feature. Correspondingly, the baryon 
mass scale of the diffusion damping is defined by Md = A^k p\)(tx jl^^^)"^ l?>^ which reduces to 

Md = 1.4 X 10^^(1 + z)-9/2(i _ yp/2)-3/2(l - /,)3/4(l^^/i2)-i/2 (20) 



- 8 - 



For the perturbations with the baryon mass scale M, this damping occurs at redshift 

1 + ZD = 1.1 X 106(1 - /,)V6(1 _ yp)-V3(f^^/,2)-l/9(^/^^)-2/9 (21) 

Here we have assumed this regime to be sufficiently before recombination. 

Now let us examine the validity of the tight coupling approximation for the baryon-photon 
system. When the photon mean free path crosses over the wavelength of the perturbations, the 
tight coupling approximation for the baryon-photon system cannot be applicable anymore. We 
define this breaking scale of the tight coupling approximation by l/k^'^ = X^^^/2tt = l/neOT- 
The corresponding baryon mass scale is defined by Mbr = 47r/?b(7r/A;g^'^)^/3, which reduces to 

Mbr = 2.9 x + z)-\l - y^/2)-\VL^h')~^ Mq, (22) 

where we have used = 1 assuming the regime sufficiently before recombination. Accordingly, 
the tight coupling approximation breaks down at redshift 

1 + ZBR = 3.8 X 10^(1 - yp)-i/2(17b/i')-'/'(M/M0)-i/6 , (23) 

for the perturbations with the baryon mass scale M. Extrapolating equation ( p2D until the decou- 
pling epoch {z ~ 1000), we obtain Mbr = 2.9 x 10^ ((1 + z)/1000)"^ (1 - yp/2)-2(17b/i^)"^ Mq. 

In Figure |2| we show the physical baryon mass scales described in the above. As can be read 
from this figure, the tight coupling approximation breaks down after the diffusion damping but 
sufficiently before recombination on small scales. Next we study the effect of this tight coupling 
breaking on the evolution of density fluctuations. 



3.3. Evolution of baryon fluctuations by the terminal velocity 

In Figure [l|, which is obtained by numerical computations, we flnd the power-law growing 
phase of baryon fluctuations after the diffusion damping but sufficiently before the decoupling on 
small scales. This evolution of baryon fluctuations is characterized as follows. Equation (12) has 
the formal solution 

aVbiv) = /"^^'"(^Qi + ^Cfc'^^' + k^)e-^''^'^''''\ (24) 

where Trf (77', ry) = J^, dfjf{fi)/ R{fj). Since the diffusion damping erases baryon-photon perturbations, 
the dominant term in the integrant of the right hand side of equation ( p^ is the gravitational force 
term. Therefore we can approximate aVb{rj) ~ Jq d-q'ak^e^'^''^'^' Before the decoupling, the 
optical depth T(i{'r]',r]) is a quite large number, and the time variation of e~'^''^'^' '^'^ is extreamly 
rapid. Thus we approximate as e""^"^ ~ g-iv-v'^M/Riv) ^ and neglect the time dependence of 
5', which derives 

Vbiv) ^ . (25) 

r 
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Fig. 2. — Physical mass scales of baryon perturbations for several different cosmological parameters. In this figure, 
'Horizon' (short dashed lines), 'Jeans' (solid lines), 'Silk' (dot-dashed lines) and 'tight couple' (dotted lines) represent 
the horizon scale, the Jeans scale, the diffusion (Silk) damping scale, and the breaking scale of the tight coupling 
approximation of baryon and photon fluids. Long dashed lines are the lines on which the Compton energy time is 
equal to the sound oscillation time scale, i.e., rjEkCc = 1. 



We can verify that this relation holds quite well during the power-law growing phase as is shown 
in Figure ^. We will discuss this figure in detail later in this subsection. The above relation is also 
obtained by equating the gravitational force due to the potential of CDM fluctuations, /c^', and 
the friction force due to the interaction with background photons, fV\^/R, in the right hand side of 



equation (12). Thus this relation implies the balance of these two forces for baryon perturbations. 
Therefore this baryon velocity can be referred as the terminal velocity. This is the result of the 
breaking of the tight coupling approximation. After the epoch ^ 1, the baryon velocity is induced 
by the gravitational free fall, i.e., Vf,(7?) = (^J^^ d-q' ak"^ + adVh{r]d)^ . Therefore the power-law 
growth by the terminal velocity lasts by the end of the baryon drag epoch. 

Let us discuss this growing feature of the evolution by the terminal velocity more quantitatively. 
In the matter dominant stage, the time variation of the curvature perturbation can be neglected. 
In this case, baryon density fluctuations are given by dh{r]) — —k'^'^{7]) p drj'R/f, from equations 
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Fig. 3. — Evolution of baryon density fluctuations and the visibility function around the decoupling epoch. The 
upper panel shows the evolution of baryon density fluctuations with the baryon mass scale 10® Mq which corresponds 
to fc = 96.5Mpc~^ (see panel (c) of Figure 1). The cosmological parameters are same as Figure 1. In this panel, 
'num' (a solid line), 'term V (a dashed line), and 'free fall' (a dotted line) represent fully numerical computation, 
the semi-analytic computation with setting Vt = k^R/T, i.e., fluctuations induced by the terminal velocity, and the 
semi-analytic computation with setting Vh to evolve via gravitational free fall after the decoupling. The lower panel 
shows the corresponding visibility functions. Note that the last scattering epoch, which is the location of the peak of 
the visibility function, of photons (a dashed line) and baryons (a solid line) are not completely the same. We refer 
the latter as a decoupling epoch in this paper. It is shown in this flgure that the semi-analytic calculation starts to 
deviate from the numerical calculations at which the decoupling process takes place (see the thickness of the visibility 
function) . 



(|Tl|) and (H). Assuming Xg = 1, we have R/f = 3.0 x lO'^Ob^i /(^e (Tt)- Then we can derive 

looo y/y ^ 

It should be noticed that this relation is applicable between the breaking epoch of the tight coupling 
and the decoupling epoch. The ratio of baryon density fluctuations to CDM ones is written as 



Kiri) ^ -0.32(1 - /„)(! - !/p/2)-'(!i„ft2)3/2 ( I^V" ( jL]\nn) . (26) 



|M.o.7xlO-(!l„-«,H-V'.(i!5!^)"\ (27) 
where we took yp = 0.23, since CDM density fluctuations can be written as 

-5.4 X 10(1 -«^(i^)(A)%(,). (28) 
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which is obtained from the relation k^"^ ~ — 47rG(a/ao)^/Oc5c- From equation (^), it is found that 
the ratio of baryon density fluctuations to CDM ones does not depend on the wave number k before 
the decoupling epoch. 

Since the above estimation is rather crude, we investigate the behavior around the decoupling 
epoch in more detail. Equation ( |2^ ) is a good approximation for 6h before the decoupling epoch. 
However, we cannot employ equation ( p6D during the decoupling epoch, because the decoupling 
process occurs in the finite time duration and our assumption of the fully ionization breaks down. 
Figure ^ shows the evolution of baryon fluctuations and the corresponding baryon and photon 
visibility functions (HS96). In the upper panel, the solid line shows the evolution obtained from 
the fully numerical computation. The dashed line shows the evolution of the baryon fluctuations 
with setting the baryon velocity to be the terminal velocity, i.e., Vf, = k^R/f. The dotted line 
shows the evolution with setting that Vf, = k^R/f until the decoupling Zd defined in HS96, and 
setting Vb to follow the free fall solution, i.e., Vb = —Vb{a/a) + A;^' after the decoupling. From this 
figure, it is apparent that our simple picture such that the small scale baryon fluctuations grow via 
the terminal velocity until the decoupling epoch, and they grows via gravitational free fall after 
that, works remarkably well before the decoupling epoch and works fairly well after the decoupling 
epoch. However, there appears a small discrepancy during the decoupling epoch, i.e., within the 
thickness of the visibility function (see the lower panel). The evolution of baryon fluctuations 
accelerates as the decoupling process begins. However the actual acceleration is not as fast as 
the one which is given by the terminal velocity (see the dashed line). This is because the friction 
between baryons and photons is getting smaller as the universe becomes transparent. If the friction 
force and gravitational force were still balanced, unrealistically large velocity would be required. 
Therefore we overestimate the velocity if we employ the terminal velocity during and after the 
decoupling epoch. This is why the dashed line exceeds the numerical one (the solid line) at the 
decoupling time which is the peak location of the baryon visibility function. Moreover, we have 
entirely ignored the friction (or drag) term after z^ in our semi-analytic calculation. Therefore 
we expect the dotted line is larger than the numerical one as is shown in Figure ^. After the 
decoupling epoch, however, the dotted line gradually overlaps with the numerical one. We note 
that this overlapping occurs earlier if we include the drag effect by photons after z^ in our free fall 
calculation. 

On the other hand, if we employ equation (p6|) by the decoupling time z^, we underestimate 
the amplitude of baryon fluctuations. Therefore the actual ratio between baryon and CDM density 
fluctuations is lower than the value of equation (p7|). Nevertheless, baryon density fluctuations have 
the amplitude of 0(10^^) ^ 0(10^^) relative to CDM fluctuations on small scales at the decoupling 
epoch. As we have already mentioned that this ratio does not depend on the wavenumber k. These 
features can be seen in Figure [l|. Since this ratio is fairly small, the assumption by HS96, that is to 
set the baryon density fluctuations to be zero after the damping epoch, is verified. Therefore we can 
employ their formulas of the evolution of CDM density fluctuations although some modifications 
are required on very small scales where the wave length of baryon perturbations is smaller than the 
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Jeans scale after the decoupling epoch. 



4. Evolution of baryon density fluctuations after the decoupling 

In this section we investigate the evolution of density fluctuations after the decoupling epoch. 
Baryon density fluctuations are catching up with CDM density fluctuations after the decoupling, 
as is shown in Figure |^. The important physical scale which is relevant to this feature is the Jeans 
scale. We define the Jeans wavelength (wave number) after the decoupling by Aj'^'^'^ = 27: /k^^^ = 
yjTiCs^/Gpm. Then we have /cj = 1.3 x 10^(1 + zY/'^{-iiT^)-^/'^{VtQh'^Y/^ Mpc^\ where Tb is the 
baryon matter temperature in unit of Kelvin, and ^ = 1 — 37/p/4. Correspondingly, the baryon mass 
scale is defined by Mj = 47ryO;,(7r/A;j'^^'')^/3, which is written as 

Mj = 1.5 X 10^(1 + z)-3/2(^^7|^)3/2(j^^;j2)-3/2^^^2 _ (29) 

Thus the epoch when the Jeans scale becomes smaller than the scale of fluctuations with M is 

1 + t = 6.1 X lQ'^-fiT^{noh^)~^{VL\,h^ f/^ {M/Mq)-^'^ . (30) 

As is shown in Figure |2|, the Jeans scale after the decoupling has the plateau. In this stage, the 
energy transfer between background photons and the baryon fluid is effective through the residual 
electrons, and the baryon temperature follows the photon temperature. As the universe expands, 
however, the energy transfer time rises above the Hubble expansion time. After that epoch, the 
baryon matter temperature cools adiabatically and drops as Tb tx a~^, which derives the decrease 
of the Jeans scale. Thus the Hubble expansion time scale and the Compton energy transfer time 
scale is equal at the broken corner of the plateau of the Jeans scale. This epoch is estimated as 
(1 + z) ~ 1000(S7b^^)^''^ ( [Peebles 1993 ). It is interesting that this redshift depends only on the 



parameter ilb^"- The line of "q^kce = 1, on which the Compton energy time is equal to the sound 
oscillation time scale, crosses at the broken corner of the plateau of the Jeans scale. Thus this cross 
point of two lines is the point when the sound oscillation time, the Hubble expansion time and the 
energy transfer time become all the same ( [YSS| ). 

This constant value of the Jeans scale, which is the maximum Jeans scale at the decoupling 
epoch, gives a characteristic scale of baryon perturbations. From the definition, this Jeans scale is 
written as /cjp = 9.0 x 10^(i7o^^)"'^^^Mpc^"'^. Correspondingly, the baryon mass is 

Mjp = 5.0 X lO^(17o/i^)^^/^Ob/i^M0 . (31) 

Here we have set = 0.23 and = 0.405, which we use in the hereafter. According to the usual 
picture of the Jeans oscillation, baryon perturbations whose scale are smaller than this characteristic 
scale are kept from growing by the baryon thermal pressure. On the other hand, the density 
fluctuations larger than this characteristic scale can grow by the gravitational infall into CDM 
potential wells, and catch up to the CDM density fluctuations. The time scale of this infall process 
is described by the free fall time. 
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By using the fitting formula, Tb = 4.5 x 10~^ {1 + z)"^ {i}\^h'^)~'^^^ K, which represents the baryon 
matter temperature after the energy transfer between baryons and photons becomes ineffective 
(YSS), we find that the epoch when the perturbations with the wave number k (> fcjp) cross the 



Jeans scale is at redshift 

1 + = 2.9 X 10^(A; Mpcy^noh^fl^h^f/^ (32) 

This epoch is rewritten in terms of the scale factor as 

* = afVoeq = 8.2 x IQ-^k Mpcfin^h^y^^^ (33) 

This is the second Jeans crossing time for the perturbations with k > kjp. These perturbations 
once cross the Jeans scale at zj of equation ([T^ ) before the decoupling epoch and starts to oscillate 
as acoustic waves. Then again they cross the Jeans scale at Zj"* and are catching up with the CDM 
density fluctuations due to gravitational infall. 



5. Transfer Function 

It will be convenient to recast the evolutionary effects in terms of a transfer function of density 
fluctuations. In the present paper, we consider the transfer function of matter density fluctuations 
which are defined by 

5m = + , (34) 

with = r^o ~ f^b- The transfer function can be defined by 

TUa,k) = .. , (35) 

where 5t is the total transfer function with radiation contributions. In the matter dominate stage, 
6m — 6t- Using the large scale solution (HS95), we obtain 

l,mM<..<^)= (l + ^/.)(l + ^/.)"5(£-)'*(0,.)AM, (36) 

where Di{a) = 2/3 + a/oeq. In general, the transfer function is a function of time. However, once 
the pressure term of baryon perturbations becomes negligible at some scale, 5^ oc Di{a) (Peebles 



1980; see also equation (|A8D ). Then the transfer function (35) becomes independent of time at this 
scale. Therefore we can assume the time invariance of the transfer function after the decoupling 
epoch on scales larger than the maximum Jeans scale after the decoupling, i.e., k < kjp or after 
Zj*"* on k > fcjp. 

While the transfer function by Bardeen, et al. (1986, hereafter BBKS) is familiar, the effect 
of the baryon fraction is not taken into account. This transfer function is not applicable for 
models with non-negligible baryon fraction at A; ;^ 0.1ft\^^^'^ftQ^'^h'^Mpc~^ . The effects of the baryon 
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Fig. 4. — Transfer functions Tm{k) at present. They are obtained by Numerical calculations (solid lines), analytic 
fitting (dashed lines) and BBKS fitting (dotted lines), respectively. 

fraction on the large-scale power spectrum have been investigated by Holtzman (1989), Peacock & 
Dodds (1995), and Sugiyama (1995). The effects of the baryon fraction on the transfer function is 
incorporated by simply scaling of so-called F factor of the BBKS transfer function However, errors 
arise on scales smaller than the galaxy scale. In Figure the fitting formula of the transfer function 
and the numerical computation described below are compared with the BBKS transfer function. 
Here the deviation is exaggerated by considering the cosmological models with high baryon density 
fraction. HS96 have investigated the CDM transfer functions on very small scales. Their result 
is summarized in the first part of Appendix B. In their paper, the effect of the baryon fraction is 
taken into account and they have given the transfer function in an analytic form. Though it is 
rather complicated equation, it has a high accuracy even for high baryon density models. However, 
there are following limitations in their transfer function. First, their transfer function is limited 
only on very small scales {k ^ IMpc), and it is not applicable at the scales larger than the galaxy 



*In a recent paper by Eisenstein & Hu (1997), they obtained better scaling, which includes k dependence, to fit 
intermediate scale (~ lOOMpc). Moreover, they presented very detailed and precise, but complicated, fitting formula 
of the transfer function on large scales {k <^ lMpc~^) with taking into account the acoustic oscillations of the photon- 
baryon fiuid before the decoupling epoch. Although their paper appears after writing up this paper, we decide to 
mention their transfer function in the following sessions. 
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scale. Secondly, they do not take into account the Jeans oscillation of baryon fluctuations after 
the decoupling epoch. Therefore their transfer function cannot apply at the wave number larger 
than kjp. Here we propose an analytic transfer function which is useful for the entire region of 
wavelength. 

The transfer function which we propose is following: 

k) = ln(l + 2.34j(l + Kj)/(l+j)) 

2Mq[l + 3.89g + (16.1g)2 + (5.460^)3 + (6.71g)4]^/^ 



with 



q = 0.951(al7e/f^o)"'/'^e|7 , (38) 

K = om9{anjno)^^^i3, (39) 



where 62.7 = ro/2.7K, a and /3 are defined in equation (B2), and f^{a,k) in equation ( p7| ) in 
Appendix B. Except for f^{a, k), the transfer function is obtained by incorporating the small scale 
transfer function in HS96 into the familiar formula of BBKS with scaling. The large scale behavior 
of the transfer function at ~ O.lMpc"^ is successfully expressed by BBKS transfer function with 
simple scaling ( ^ugiyama 1995| ). However, the small scale behavior is not satisfactory. It is required 
that the transfer function approaches to the form in HS96 in the small scale limit. The term in the 
logarithm in equation (|37|) plays a role to change from large scale formula to the small scale one. 



Let us briefly comment on the correction factor f^{a,k), which is described in detail in 
Appendix B. As discussed in §4, the maximum Jeans scale after the decoupling epoch is ex- 
pressed by /cjp. The baryon thermal pressure keeps baryon fluctuations from growing even af- 
ter the decoupling epoch on the scales k > kjp. When baryon fluctuations can not grow and 
remain the small amplitude, the growth rate of CDM density fluctuations is suppressed. In 
this stage, the growth rate of CDM fluctuations is roughly in proportion to (1 + a/oeq)""^ with 
ai = (1 — y^l + 24(r2o — r2b)/f^o)/4. Eventually this effect decreases the transfer function on sales 
smaller than A;jp after the decoupling. In order to describe this effect, we introduce the window 
function /m(a, k). 

In Figure 5, the fitting formula of the transfer function is compared with the numerical com- 
putation, normalized to the BBKS transfer function. Our fitting formula reproduces the numerical 
result quite well for the scale, IMpc^^ ^ ^ ^ lOOMpc^^, within a few percent for 0.5 <^ h <^ 0.8 
and Q\^/Qq <^ 0.5. For the smaller scales k ^ lOOMpc"^, where fm{<i-,k) becomes important, our 
fitting function works within ~ 10% accuracy at z ^ 100. For the larger scales k <^ O.lMpc"^, the 
numerical transfer function shows the feature with bumps and wriggles due to the baryon-photon 
acoustic oscillation before the decoupling. The larger the baryon fraction is, the bigger this effect 
becomes. As we do not take this effect into account, the difference between the fitting formula and 
the numerical result becomes large. However, the fitting formula still designs to cross the center 
of these oscillations. Therefore we expect fairly good fit with the observational quantities, such as 
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(78. The accuracy of the fitting formula is checked in the next section in more detail by computing 
statistical quantities. 
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6. Simple Demonstrations 

In this section, we first show the accuracy of the fitting formula by calculating the following 
statistical quantities of the large scale structure the universe. Here we consider the CDM cosmology 
with the Harrison-Zeldovich initial density power spectrum. For low density universe models, we 
introduce the cosmological constant to keep the flat geometry. Let us introduce which is defined 
as 

1 f°° 

a\{a) = — dkk'^P{a,k)W{kRf , (40) 
zvr^ Jo 

where P{a,k) = {5m{a,k)'^) is the power spectrum of the linear density perturbations at time a, 
and W{x) = 3(sinx — xcosx)/x'^ is the top-hat window function. To determine the amplitude of 
the power spectrum, we use the COBE normalization by Bunn &: White (1996). According to their 
results, we can rewrite equation (|^) as 

„ ^ 6l — T,„ a.k?W kR?. 



= I «;Mj^«i„ -(^-j T„(„..)^HW. (41) 

with 6h = 1.94 X io-5j7o'°-'^^^"°-°^^°^o for A-models. We note that Di{a) in equation (|l|) represents 
the growing mode solution for matter perturbations taking A-term into account, which differs from 



D,{a{l + z)) = -^Jno{l + z)3 + 1 - / da'i——-^—-—] , (42) 



(A8) at z ^ a few due to the A-term. For A-model, Di{a) is written as (see e.g., Peebles 1980) 

a' 

\no + a'^{l-a 

which is normalized as Di{a) ~ 1/(1 + z) at z ^ 1. 

The observational quantity ag is defined as erg = (^R=8h-^Mpci^ = 0). In Table ||, we show the 
value of CTg by numerical calculations comparing with various transfer functions. In this table, , 
(Tg, (Tg , (T^^^^ and o"g^ denote values obtained by the numerical calculation, our fitting formula, 
the empirical scaling by Sugiyama (1995), the original fitting formula by BBKS and the scaling 
by Eisenstein & Hu (1997), respectively. 

It is shown that our fitting formula together with that obtained by the empirical scalings by 
the shape parameter F by ( Sugiyama 1995| ) and by Eisenstein & Hu (1997), well reproduce the 



numerical result. This is because that the fitting formulas are designed to cross the center of 
oscillations. 

In order to check the accuracy on small scales, we have calculated the quantities, CR=ih-iMpc('2^ = 
3), (yR=o.ih-mpc{z = 5), CTR=iQ-2h-iMpc{z = 10), and (Ti?=io-3h-iMpc(^; = 10), for various cosmo- 
logical models. The results are shown in Tables From these results, we can find that our fitting 
formula works very well on small scales. 



^ Here we used Tcfi{k) in their paper in stead of the usual F- factor (= floh). 
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Fig. 6. — Mass fractions of the bound objects F{> Mr). The cosmological parameters taken in this figure is Qo = 0.3 
and h — 0.6. In each panel the baryon fraction is different, which are: (a)f2b ~ 0; (b)S7b — 0.05; (c)r2b = 0.1; and 
(d)!:!b = 0.15. 



Next we briefly demonstrate a cosmological implication of our results. The study of the early 
formation of collapsed objects and the thermal history of the high-z universe is one of the most 



important issue in cosmology and galaxy formation, (e.g., Fukugita &: Kawasaki 1994 ; Gnedin fc 
Ostriker 1997 ; Ostriker &: Gnedin 199£| ; Haiman & Loeb 1997). The statistical arguments allow 
us to investigate such early history of the formation of small scale cosmic objects in the high-z 
universe without big numerical calculations. The Press-Schechter theory (Press & Schechter 1974) 
is among the most simple theory to discuss the statistics of the gravitationally collapsing objects. 
According to the Press-Schechter formula, we calculate the fractions of the gravitationally bound 
system, and show the effects of the baryon fraction on the formation rate of the nonlinear objects 
in the high-z universe. 

The Press-Schechter theory predicts that the mass fraction which is gravitationally bounded 



objects above a given mass M/j at cosmic time a becomes (e.g., Padmanabhan 1993 ), 



F(> Mr) =erfc(- 



(43) 



where a{Mji, a) 
matter energy density, and 6, 



\V2a{MR,a) 

crR{a) is define by equation (|4l|), Mr = AirpoR^/S, po is the spatially averaged 



1.69. Figure g shows the mass fraction F{> Mr) calculated from 
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equation (|43|), for the cosmofogical model JIq = 0.3, h = 0.6 with different baryon density fractions, 
(a) is the case that the baryon fraction is neglected, (b) is the case with including the baryon 
fraction = 0.05. (c) is the case = 0.1, and (d) is = 0.15. It is apparent that the formation 
rate of cosmic objects is delayed in a universe with large baryon fraction. 

7. Summary and Discussions 

In this paper, we have carefully re-investigated the evolution of the small scale density fluc- 
tuations in an expanding universe, motivated from the studies of cosmic structure formation in 
the high-z universe. Under the assumption of the hierarchical clustering scenario, the structure 
formation takes place from smaller scales. Therefore such careful investigation of the evolution of 
density fluctuations on small scales is necessary. In the universe with non-negligible baryon frac- 
tion, the evolution of density fluctuations is very complicated because baryon perturbations affect 
CDM perturbations through the gravitational interaction. 

Interestingly, baryon density fluctuations show the varieties of physical processes caused by 
the interaction with photons which are summarized in Figure ||. Discussing the physical scales 
characterizing these processes, we have studied the evolution of the small scale baryon density 
perturbations. In this investigation, we have found the growth of baryon density fluctuations on 
small scales after the diffusion damping and before the decoupling epoch, which is related with 
the breaking down of the tight coupling approximation between baryons and photons. We have 
referred this growing mode as the terminal velocity mode. 

We have re-analyzed the total matter transfer function. In the universe with non-negligible 
baryon fraction, the transfer function is changed from the familiar transfer function by BBKS, 
for k ^ O.mo^/^no^/'^h'^. Extending the previous fitting formulas by Sugiyama (1995) and IIS96 
with including the Jeans effect after the decoupling epoch which causes the damping on very small 
scales, we have presented the analytic transfer function in a way of modifying the BBKS transfer 
function. This fitting function is rather complicated, but its accuracy is very high particularly 
around the galaxy scale IMpc^^ ^ k <^ lOOMpc^^. This fitting function is applicable for the 
entire range. Even for k ^ lOOMpc^^ it works within ~ 10% error. Unless one would like to 
trace the oscillation feature on intermediate scales 0.01Mpc~^ ^ ^ ^ O.lMpc"^ which is caused by 
the acoustic oscillation by photon-baryon fluid before the decoupling epoch (in this case, one need 
to employ pure numerical calculations or the newly suggested fitting formula by Eisenstein & Hu 
(1997)), our fitting works fairly well for the entire range of k ^ lO^Mpc"^. 

When discussing the formation epoch of small cosmic objects, the amplitude of density fluctu- 
ations is the important factor. In the universe with large baryon density fraction, the amplitude of 
total matter perturbations is decreased on small scales. Then the formation of the cosmic objects is 
delayed. Moreover, the evolution of total density fluctuations on these scales does not quite follow 
the usual growth rate Di. We have described this time dependence of the transfer function as /m. 
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In a low density universe, the baryon density fraction becomes large relative to the CDM density 
fraction. Then this effect becomes significant, as demonstrated in §6. 

Finally, if once the early re-ionization occurs, baryon density perturbations evolve in a different 
way. This effect is beyond the scope of this paper. But a recent paper by Chiba, Kawasaki & 
Sugiyama (1997) showed that the baryon and total matter evolution are affected by the early 
re-ionization in case of high baryon fraction models. 
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A. Evolution of Matter Perturbations 

In this Appendix, we briefly review the analytic formulation of the evolution of the matter 
perturbations on small scales obtained by HS96. Let us consider the evolution after the diffusion 
damping. If we set photon and baryon fluctuations to be negligible small (this assumption is valid 
by the argument in §3.3 even though there exists small fraction of baryon fluctuations after the 
damping epoch due to the terminal velocity), only CDM contributes to the gravitational potential. 
Which permits the equation of CDM perturbations to be written as 

£6, 2 + 3y d5, ^ 3 . 
dy^ ^ 2y{l+y) dy 2y{l + y) ^ ^ 

where y = a/ocq and f^c = — J^b- The solutions are given in terms of the Hyper geometric 
function 

U^ = {l + y)-°»2i^i (tti, ai + 1/2, 2a, + 1/2; 1/(1 + y)) , (A2) 



where i = 1, 2 and = (1 it y/1 + 24r2c/f^o)/4 with — and + for i = 1 and 2, respectively. 

With considering the matching conditions to the solution outside the horizon in the radiation 
dominated stage, the solution for 5c after the horizon crossing is obtained. 



5ciri,k)=lM0,k)[AiUi{r])+A2U2iv)), (A3) 

where 



^ r(ai)r(ai + i/2) n^^ 

' r(2Qi + 1/2)(V'(2q2) - V(2ai)) V2 



-'2 



+ ^^(1) - V'(a2) +ln2^, (A4) 



OH- 

A2 is obtained by replacing the subscripts 1 2, /i and I2 are obtained from the numerical fitting 
formulas 

h = 9.11(1 + 0.128/^ + 0.029/^2)^ (A5) 
h = 0.594(1 - 0.631/^ + 0.284/^^), (A6) 

QH ^ l + 7l + 8(/c/fceq)^ 

and 'ip{x) is the poli-gamma function. 

After the epoch that the baryon pressure becomes negligible, i.e., the drag epoch or the 
second Jeans crossing epoch Zj"* depends on scales, baryon fluctuations evolve through CDM 
gravitational potential well. In this stage, it is well known that the matter fluctuations (5m, defined 
by equation (p^, follow the growing and decaying solutions (Peebles 1980) 

Di = ^ + y, (A8) 



D2 = ^(2 + 3y)ln 



:i + y)V2 + i 



(l + y)V2_l 



(A9) 
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These solutions correspond to Ui and U2 with ai\n^=ng = (1 it 5)/4 . Matching solutions at 
y = Uc = Vd or = yj"*, we obtain for the growing mode solution as 



(AlO) 



y=yc 



with 



Gi 



Do 



D1D2 - D1D2 



Go 



Do 



D1D2 - D1D2 



(All) 



B. Effect of the baryon thermal pressure on the Transfer Function 

In this Appendix we describe some details which are important to construct an analytic transfer 
function on very small scales, k > IMpc^^. We first note the small scale transfer function obtained 
in HS96. According to their result, the fitting formula of the CDM transfer function is written as, 

ln[1.8/3g] 

Tc{a,k) = a-^j^^, (Bl) 



with 



/3-i = l + 6i((J^cM)''-l), (B2) 



ai = {Aemoh^f^^il + (32.ir?o/i^)"°-^^^) 
a2 = {l2moh^f^^\l + (45.0^0/^')^°-''') 
bl = 0.944(1 + (458.00o/i^)"°-^°^)"\ 
62 = (0.395^^o/l^)"°•°^^^ 



and q = A;(Mpc i^^h ) (To/2.7K) . Though this formula is rather complicated, it works to 1% 
accuracy. We should notice that this formula is applicable only on small scales, k ^ IMpc^^. 







Let us describe the function /m(a, A;) in equation (|37D. As we have shown in Appendix A, 
6ja grows in proportion to Di(a) after the pressure becomes ineffective. On scales larger than the 
Jeans scale at the decoupling epoch, i.e., k < /cjp, the photon pressure becomes ineffective after the 
decoupling epoch. Then 6m evolves in proportion to Di{a). Once the matter fluctuations become 
to follow this growing mode solution, the transfer function becomes independent of time. Therefore 
on these scales, the transfer function does not change its shape after the drag epoch. On scales 
k ^ /cjp, on the other hand, the baryon thermal pressure keeps the perturbations from growing 
even after the decoupling of baryons from photons. As we have seen in Appendix A, the growth 
rate of CDM fluctuations is suppressed in the situation that there exists a homogeneous matter 
(baryon) component besides the CDM component. This effect alters the transfer function on small 
scales k ^ kjp even after the decoupling epoch. The function fm{a,k) is multiplied in order to 
describe this effect. 
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In order to find the analytic form of fmia, k), we rewrite 

T^{a, k) = rm(ad, k) 3"/"' = T^ia^, k)f{a, k), (B3) 
-/m(ad, k) 

where Od denotes the scale factor at the decoupling epoch. Thus /(a, k) is defined as the ratio of 
the transfer function at a{rj) to that at the decoupling time. From the definition of the transfer 
function (|35[), we can rewrite 

f{a,k) = — — . (B4) 

Di{a) dia{ad,k) 

We first consider the case k < kjp. This case is trivial, because 5^ grows in proportion to 
Di{a) after the decoupling, as described before. Then we get f{a,k) = 1. 

We next consider the case k > kjp. These perturbations stay inside the Jeans scale until 
aj'^^/fleq = Uj'^^i which is defined in equation (33). Assuming that baryon density fluctuations do 
not contribute to the gravitational potential in this regime, we approximate the growth rate of 
the CDM density perturbations by the growing mode solution Ui{a) in equation ( [A2D . After the 
decoupling epoch, we can assume the matter domination of the universe, i.e., a/oeq ^ 1. Hence 
we further approximate Ui{a) ~ (1 + a/aeq)"^ with ai = (1 — ^1 + 24r2c/f^o)/4 as is shown in 
Appendix A. Then we can write 

where y = a/ueq and = a^/a as are defined before. After the Jeans crossing a/oeq > yj"*, we 
approximate that (5m follows the growing mode Di{a) and that the solution is given by equation 
( AlOP with replacing yc by yj"*. Then we can write 



f{a,k) ^ ^h^(G,{a)Ui{a)-G2{a)Ui{a)) 



a/aGq=J/j" 



3-2ai |^l+ad/aeqy+"i 

' ' , (a/aeq > yj ). (B6) 



5 V 1 + 

Here we have assumed that the amplitude of 5^ is small at yj"*. 

Summarizing the results, /(a, k) is evaluated 

1 {k < kjp), 

f{a,k)^{ ((l + yd)/(l + y))"'^^ (/c> /cjP and y < y;?'^*), 

((1 + yd)/(l + yr))°''''(3 - 2ai)/5 {k > A;jp and y > yf 



In the case when the baryon fraction is negligible, i.e., Ob *C ^O) we have ai = —1 and /(a, k) = 1. 
On the other hand, in the cosmological models with the non-negligible baryon fraction, the growth 
rate of CDM fluctuations is suppressed even after the decoupling. And the matter transfer function 
is changed. 
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Combining above /(a, A;)'s on difference scales, we give tfie next fitting function, 



fm{a,k) 



+ 



(fcAjp)l-3 + 1 (A;jp/fc)l-3+0.9 

■i 



l + Vd 

(y/yOut)1.3 + i\^ 1+y 



ai+1 



3-2qi / 1 + yd 



with 



(yOUt/y)1.3 + 1 5 yi+yO 

fejp = 350(1 + 50 (yd/y))(^^o/i')'/'Mpc-\ 



..out 



ai+1 



, (B7) 



(B8) 
(B9) 



and y = a/ogq = 2.4x 10^(l+2;)^^no^^. Here we determine numerical constant factors to reproduce 
the numerical computations. The reason why we employ /cjp instead of fcjp is following. Even on 
scales smaller than /cjp, the thermal pressure cannot entirely prevent baryon fluctuations from the 
growth because of the existence of the CDM potential. A certain amount of the growth rate has 
to be taken into account at ;^ kjp even though the rate is less than Di. 
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Table 1. crfl=8/i-iMpc for various cosmological models. 



no 


fib 


h 


N 


F 

"Sh-lMpc 


_s 

"sa-lMpc 


BBKS 
"sh-iMpc 


EH 
"Sh-lMpc 


0.3 


0.05 


0.6 


0.70 


0.71 


0.70 


0.95 


0.71 


0.3 


0.1 


0.6 


0.51 


0.51 


0.51 


0.95 


0.50 


0.3 


0.15 


0.6 


0.35 


0.33 


0.37 


0.95 


0.34 


0.3 


0.2 


0.6 


0.23 


0.19 


0.26 


0.95 


0.22 


0.3 


0.05 


0.8 


0.98 


1.00 


0.98 


1.36 


0.99 


0.3 


0.1 


0.8 


0.68 


0.69 


0.70 


1.36 


0.68 


0.3 


0.15 


0.8 


0.46 


0.44 


0.49 


1.36 


0.44 


0.3 


0.2 


0.8 


0.29 


0.24 


0.34 


1.36 


0.28 


1.0 


0.05 


0.5 


1.18 


1.22 


1.18 


1.31 


1.21 


1.0 


0.1 


0.5 


1.08 


1.12 


1.06 


1.31 


1.11 


1.0 


0.5 


0.5 


0.50 


0.46 


0.41 


1.31 


0.48 


1.0 


0.05 


0.8 


1.83 


1.92 


1.87 


2.07 


1.90 


1.0 


0.1 


0.8 


1.65 


1.76 


1.68 


2.07 


1.73 


1.0 


0.5 


0.8 


0.74 


0.66 


0.62 


2.07 


0.70 



'^a^ ,a^^^^ , and cr™ are calculated by the full numerical computation, 

using our analytic fitting, an analytic fitting by Sugiyama (1995), the one by BBKS 
and the one by Eisenstein & Hu (1997) with employing their empirical scalings of 
the shape parameter F, respectively. 
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Table 2. aj^=i,j-iMpc at z = 3 





fib 


h 


N 

°lh-lMpc 


'^Ih-lMpc 


'^Ih-lMpc 




"-"^Ih-lMpc 


0.3 


0.05 


0.6 


0.66 


0.66 


0.64 


0.92 


0.65 


0.3 


0.1 


0.6 


0.44 


0.44 


0.44 


0.92 


0.43 


0.3 


0.15 


0.6 


0.28 


0.27 


0.30 


0.92 


0.26 


0.3 


0.2 


0.6 


0.16 


0.15 


0.20 


0.92 


0.14 


0.3 


0.05 


0.8 


1.00 


0.99 


0.96 


1.45 


0.97 


0.3 


0.1 


0.8 


0.64 


0.63 


0.64 


1.45 


0.61 


0.3 


0.15 


0.8 


0.38 


0.37 


0.42 


1.45 


0.36 


0.3 


0.2 


0.8 


0.20 


0.19 


0.27 


1.45 


0.19 


1.0 


0.05 


0.5 


1.31 


1.31 


1.25 


1.45 


1.29 


1.0 


0.1 


0.5 


1.16 


1.16 


1.08 


1.45 


1.15 


1.0 


0.5 


0.5 


0.37 


0.37 


0.31 


1.45 


0.36 


1.0 


0.05 


0.8 


2.51 


2.53 


2.44 


2.86 


2.51 


1.0 


0.1 


0.8 


2.18 


2.21 


2.07 


2.86 


2.18 


1.0 


0.5 


0.8 


0.58 


0.57 


0.53 


2.86 


0.56 



, and a are same as Table 1. 



Table 3. 0-^=0. i^-iMpc at 2; = 5 



no 


fib 


h 


N 

°0.1h-lMpc 




Ih-lMpc 


BBKS 
"o.lh-lMpc 


EH 

"o.l/i-lMpc 


0.3 


0.05 


0.6 


0.94 


0.93 


0.88 


1.32 


0.89 


0.3 


0.1 


0.6 


0.61 


0.60 


0.59 


1.32 


0.57 


0.3 


0.15 


0.6 


0.36 


0.36 


0.39 


1.32 


0.33 


0.3 


0.2 


0.6 


0.19 


0.19 


0.26 


1.32 


0.17 


0.3 


0.05 


0.8 


1.47 


1.46 


1.38 


2.16 


1.39 


0.3 


0.1 


0.8 


0.91 


0.90 


0.88 


2.16 


0.84 


0.3 


0.15 


0.8 


0.51 


0.51 


0.56 


2.16 


0.47 


0.3 


0.2 


0.8 


0.25 


0.25 


0.35 


2.16 


0.24 


1.0 


0.05 


0.5 


2.26 


2.23 


2.08 


2.47 


2.16 


1.0 


0.1 


0.5 


1.98 


1.95 


1.76 


2.47 


1.89 


1.0 


0.5 


0.5 


0.55 


0.55 


0.45 


2.47 


0.51 


1.0 


0.05 


0.8 


4.83 


4.75 


4.48 


5.41 


4.63 


1.0 


0.1 


0.8 


4.13 


4.06 


3.71 


5.41 


3.93 


1.0 


0.5 


0.8 


0.93 


0.91 


0.80 


5.41 


0.84 




0- ,cr , 


^BBKS 


, and cr^^ are 


same as Table 1. 
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Table 4. o-^=io-2/i-iMpc at 2: = 10 



^0 


fib 


h 


'^l()-2h-lMpc 




"^lO-^h-lMpc 


''^lO-2/i-lMpc 


''^lO-^h-lMpc 


0.3 


0.05 


0.6 


0.85 


0.85 


0.79 


1.20 


0.80 


0.3 


0.1 


0.6 


0.54 


0.54 


0.52 


1.20 


0.50 


0.3 


0.15 


0.6 


0.32 


0.32 


0.34 


1.20 


0.29 


0.3 


0.2 


0.6 


0.16 


0.16 


0.22 


1.20 


0.14 


0.3 


0.05 


0.8 


1.37 


1.36 


1.26 


2.02 


1.27 


0.3 


0.1 


0.8 


0.82 


0.83 


0.79 


2.02 


0.75 


0.3 


0.15 


0.8 


0.45 


0.46 


0.49 


2.02 


0.41 


0.3 


0.2 


0.8 


0.22 


0.22 


0.31 


2.02 


0.20 


1.0 


0.05 


0.5 


2.28 


2.25 


2.05 


2.46 


2.14 


1.0 


0.1 


0.5 


1.98 


1.96 


1.72 


2.46 


1.85 


1.0 


0.5 


0.5 


0.52 


0.53 


0.41 


2.46 


0.47 


1.0 


0.05 


0.8 


5.12 


5.04 


4.63 


5.65 


4.79 


1.0 


0.1 


0.8 


4.34 


4.28 


3.79 


5.65 


4.03 


1.0 


0.5 


0.8 


0.91 


0.93 


0.75 


5.65 


0.79 



, and a are same as Table 1. 



Table 5. o-^^^io-s/i-iMpc at z = 10 



''^lO-3/i-lMpc ^10-3h-lMpc '-'^lO-Sh-lMpc '^IQ-^ h.-'^Mpc "^10-3 ft- 1 Mpc 



0.3 


0.05 


0.6 


1.16 


1.15 


1.15 


1.77 


1.17 


0.3 


0.1 


0.6 


0.68 


0.68 


0.75 


1.77 


0.72 


0.3 


0.15 


0.6 


0.38 


0.38 


0.49 


1.77 


0.41 


0.3 


0.2 


0.6 


0.19 


0.19 


0.32 


1.77 


0.20 


0.3 


0.05 


0.8 


1.88 


1.86 


1.86 


3.02 


1.88 


0.3 


0.1 


0.8 


1.05 


1.04 


1.15 


3.02 


1.10 


0.3 


0.15 


0.8 


0.55 


0.55 


0.71 


3.02 


0.59 


0.3 


0.2 


0.8 


0.25 


0.26 


0.44 


3.02 


0.29 


1.0 


0.05 


0.5 


3.48 


3.47 


3.16 


3.80 


3.29 


1.0 


0.1 


0.5 


2.97 


2.94 


2.63 


3.80 


2.83 


1.0 


0.5 


0.5 


0.69 


0.70 


0.60 


3.80 


0.69 


1.0 


0.05 


0.8 


8.05 


8.00 


7.30 


8.97 


7.56 


1.0 


0.1 


0.8 


6.68 


6.62 


5.93 


8.97 


6.32 


1.0 


0.5 


0.8 


1.23 


1.25 


1.12 


8.97 


1.18 


0^ ; 


a ,a , 


^BBKS 


, and a™ 


are same as Table 1. 









